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Abstract 



We calculate the Wigner quasi-probability distribution for position and momentum, (x,p), 
for the energy eigenstates of the standard infinite well potential, using both x- and p-space 
stationary-state solutions, as well as visualizing the results. We then evaluate the time-dependent 
Wigner distribution, P]y(x,p;t), for Gaussian wave packet solutions of this system, illustrating 
both the short-term semi-classical time dependence, as well as longer-term revival and fractional 
revival behavior and the structure during the collapsed state. This tool provides an excellent way 
of demonstrating the patterns of highly correlated Schrodinger-cat-like 'mini-packets' which appear 
at fractional multiples of the exact revival time. 
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^Electronic address: madlO@psu.edu 
+ Electronic address: rick@phys.psu.edu 





2 



I. INTRODUCTION 



The solution and visualization of problems in one-dimensional quantum mechanics fo- 
cuses most often on calculations of the position-space wavefunction, ip(x,t). More occasion- 
ally, problems may be solved using, or transformed into, the momentum-space counterpart, 
4>(p,t), using the standard Fourier transforms, 

i r+oo 

i>{x i t) = -=j 4>(p,t)e +i ^ h d P (1) 

and 

i r+oo 

<j>(p,t) = -j= J ^(x,t)e-^ h dx. (2) 

Connections between the classical and quantum descriptions of model systems can then be 
made in a variety of ways. For example, the quantum mechanical expectation values, (x)t and 
(p)t can be compared to their classical analogs, x(t) and p(t) = mv(t), via Ehrenfest's prin- 
ciple, e.g., (p) t = md(x} t /dt. Quantum mechanical probability densities, Pqm(x) = |?/> n (x)| 2 
and Pqm(p) = I0n(p)| 2 , can be related to classical probability distributions, calculated using 
simple "how long does a particle spend in a given dx or dp interval?' arguments, or justified 
in more detail through the WKB approximation. 

The visualization of solutions to problems in classical mechanics through a phase-space 
description, i.e., parametric plots of p(t) versus x(t), is often helpful. It is natural to wonder 
if a quantum mechanical analog of a phase-space probability distribution, a joint P(x,p) 
probability density, is a useful construct, despite the obvious problems raised by the Heisen- 
berg uncertainty principle and its connection between x and p. 

Wigner [l| was one of the first to address this issue and introduced a quasi- or pseudo- 
probability distribution corresponding to a general quantum state, ip(x,t), defined by 

i r+oo 

P w ( XjP; t) = — / r(x+y,t)iP(x-y,t)e 2i ^ n dy. (3) 



The properties of the Wigner distribution have been discussed in a number of accessible 
reviews 



applications [1 



'Perties ol tne Wigner distribution nave been discuss 

13 - |lfl| [. a naly zed for theoretical consistency ^j] - us ed for a number of physical 



discussed in the pedagogical literature [2(j] - 2£ | , and have even made 
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Similar topics 



the occasional brief appearance in quantum mechanics textbooks 
are often discussed in the context of both the signal analysis of time- var ying spectra 
(where the two complementary variables are t and u) and quantum optics [29j . [30]. 



Of all of the familiar model systems often used as tractable pedagogical examples for 
the application of quantum mechanical techniques, many of the standard potentials have 
been analyzed in detail using the Wigner distribution, including the harmonic oscillator, 
linear potential, Morse potential, hydrogen atom (Coulomb problem), and others. The most 
familiar example of all, however, the infinite square well (ISW) problem, has received less 
attention P|, [n]], j^. The Wigner function has recently been used to help explain the 
interesting patterns of probability density, P(x,t) = \ip(x,t)\ 2 versus (x, t), when plotted 
over long time periods (namely, over the interval < t < T rev where T rev is the revival time, 
to be discussed below). Such patterns have come to be known as quantum carpets 
The authors of Ref. [40| have produced video sequences of the free evolution of an array of 
Wigner functions constructed from an antisymmetric wavefunction of a wave packet and its 
mirror wave packet in order to understand the 'weaving' of quantum carpets. While such 
graphical representations of the long-term time dependence of this familiar problem are very 
interesting, the necessary background material which would allow students (or instructors) 
to approach such problems themselves has not been stressed. 

The purpose of this paper, therefore, is to thoroughly review the derivation of the Wigner 
distribution for this very accessible problem, not only for individual energy eigenstates, but 
for time-dependent wave packet solutions. This last case is important as the Wigner function 
provides a very useful tool for the visualization of the non-trivial long-term time evolution 
of initially localized wave packets, including revival and fractional revival behaviors. (For 
a pedagogical review of this topic, see Ref. and for a recent survey of the research 



literature, see Ref. (4^ 

In the next section, we provide a short, self-contained review of many of the general 
properties of the Wigner quasi-probability distribution, followed in Section lTTTl bv a collection 
of useful results for the ISW energy eigenfunctions, in both position and momentum space, 
helpful for the evaluation of the Wigner function for the ISW. We briefly review the classical 
phase-space description of the infinite square well in Sec. IIV[ and in Section El we exhibit 
the Wigner distribution, P^\x,p), for the ISW eigenstates, illustrating their derivations in 
both x- and p-space, and comparing them to their classical counterparts. In Section IVI1 we 
extend these results to general wave packet solutions of the infinite square well, focusing on 
the structure of the Wigner function for time-dependent Gaussian wave packet solutions of 
the ISW, at a number of full and fractional revivals, as well as during the so-called collapsed 
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phase, illustrating the power of this method for the visualization of time- dependent quantum 
phenomena. Throughout, we focus on providing explicit analytic results which can be used 
for further numerical or visual investigations. 



II. REVIEW OF THE WIGNER QUASI-PROB ABILITY DISTRIBUTION 

The definition of the Wigner function, 

i r+oo 

P w ( x , p]t ) = — / r(x+y,t)ip(x-y,t)e 2i ™/ h dy, (4) 

J— oo 

using position-space wavefunctions, can be easily rewritten in momentum space, using 
Eqn. (J2J), in the equivalent form 

i r+oo 

P w (x,p;t) = — / <P*(p + q,t)^p-q,t)e- 2 ^ h dq. (5) 

™ J-oo 

The Wigner distribution is easily shown to be real, 

i r+oo 

[Pw(x,p;t)}* = — / i>{x + y,t)r(x-y,t)e- 2 ™/ n dy 



oc 
oo 



/ ^(x + y,t) ^{x - y, t) e +2l ^ h dy (6) 
= Pw(x,p;t) 

by using a simple change of variables (y = —y). This is, of course, one of the desired 
properties of a probability distribution. 

Integration of Pw(x,p; t) over one variable or the other is seen to give the correct marginal 
probability distributions for x and p separately, since 

/+oo 
P w (x,p;t)dp = \*p(x } t)\ 2 = P QM (x,t) (7) 
-oo 

/+oo 
P w (x,p;t)dx = \<P(p,t)\ 2 = P QM (p,t) (8) 
-oo 



where one uses the definition of the Dirac 5 function 



1 



+oo 



5{z) = — / e lkz dk (9) 

in Eqns. (@J) or (jSJ). This property of Pw(x,p) is clearly another necessary condition for a 
joint probability density. 
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However, one can also easily show that the Wigner distributions for two distinct quantum 

states, ip(x,t) and x{ x ,t), 

i r+oo 

P$\x,V\t) = ^ I r(x + y,t)Hx-y,t)e 2 ^ h dy (10) 



7T 



— oo 
+00 



1 r+oo 

P${x,F,t) = -r \ X *(x + z,t) X (x-z,t)e 2 ^ h dz (11) 
satisfy the relation 

/+OO r+OO r) 

dx / dpP$\x,p-t)P$(x,p;t) = — Mx)\ 2 . (12) 
-00 J -00 

So, for example, if ip and x are orthogonal states so that (ip\x) = 0, it cannot be true 
that the corresponding Wigner distributions can be everywhere non-negative, as there must 
be cancellations in the integral in Eqn. ()12j) . The fact that Pw{x,p;t) can be negative is 
easily confirmed by direct calculation, for example, with simple cases such as the harmonic 
oscillator (for any n > 1 state). This feature is, after all, hardly surprising because of the 
non-commutativity of x and p encoded in the uncertainty principle. Despite this apparent 
drawback the Wigner function is still useful for the visualization of the correlated 
position- and momentum-space behavior of quantum eigenstates and wave packets. 

A useful benchmark example of the Wigner function is for a Gaussian free-particle wave 
packet, where the time-dependent momentum- and position-space wavefunctions (for arbi- 
trary initial Xq and po) can be written in the forms 



(f> {G) (p,t) = l® e -a 2 (P-Po) 2 /2 e -ipxo/h e - ip it/2 m h ( 13 ) 



lb(r\(x t) = = p ipo(x-x )/h p -iplt/2mh p -(x-x -p t/m) 2 /2(ah) 2 (1+it/tp) Q^N 

' ^ah(l + it/t ) 



which are related by Eqns. (0) and (j2J). These solutions are, of course, well-known to be 
characterized by 

(p)t=Po, Ap t = -^= (15) 

ay 2 



(x) t = (po/m)t + x , Ax t = (ha/ V2)y/1 + (t/t ) 2 (16) 

where tn = mha 2 is the spreading time. The corresponding Wigner function is easily ob- 
tained [a], using standard Gaussian integrals, and is given by 

P$\x,p; t) = ^z e - a2{p - po)2 e -(—o-pt/m)y^ (17) 



where f3 = ha. In this case, the ultra-smooth Gaussian solution does give a positive-definite 
Pw{x,p;t) and it is known 4^ that such solutions are the only forms which give rise to 
non-negative Wigner functions. 

Another result which will prove useful in what follows is the expression for the Wigner 
function for the case of a linear combination of two 1-D Gaussians, characterized by different 
values of xq and po- For example, if we assume that at some instant of time we have 



ip( A > B \x) = ltp(G)(x;x A ,p A ) + 5iP( G )(x;x b ,Pb) 



(18) 
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e -(x-x A ) 2 /2(3 2 e ip A (x-x A )/h 
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+ 5 



e -(x-x B ) 2 /2f3 2 e ip B (x-x B )/H 
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the corresponding Wigner function is given by 



w 



where 



nh 



^2 e -* 2 (p-p A ) 2 e -(x-x A ) 2 /p 2 + ^2 e -a 2 (p-p B ) 2 e -(x-x B ) 2 /p 2 



(19) 



+ 2e -a 2 {p-p) 2 e -(x-x) 2 IH 2 Re { lS * e i(* 



A p B ~x B p A )/h -i(x A -x B )(p-p)/h i(p A -p B )x/h^ 



_ Xa+Xb 
X 



and 



_ Pa+Pb 
p = 



- ■ 2 - < 20 > 

In this case, the Wigner function is characterized by two smooth 'lumps' in phase space, 
corresponding to the values of (xa,Pa) and (xb,Pb) of the individual Gaussians, but also 
by an oscillatory term, centered at a point in phase space defined by the average of these 
values; the oscillations appear in the x variable (if pa ^ Pb), the p variable (if xa # %b) 
or both. (A similar expression for the case where pa = Pb — is given in Ref. [27i ]; see 
also Refs. 2^ and j^] for related examples involving time-frequency analysis and coherent 
states, respectively.) 

For a single stationary state or energy-eigenfunction solution, corresponding to a quan- 
tized bound state, given by 

j} n {x,t) = u n (x)e- iE ~ t " i (21) 

where u n (x) can be considered as a purely real function, the Wigner distribution is time 
independent 



P${x,p-t) 



1 



u n (x + y) u n (x - y) e M ™' n dy = P$ } (x, P ) . 
For such solutions, another change of variables argument suffices to show that 



(22) 



P%\x,^>)=P%>(x,+p) 



(23) 



which corresponds to the classical result that a particle undergoing bound state motion will 
spend equal amounts of time with +p (motion to the right) as in the — p direction (going 
to the left). A more familiar version of this is often discussed in textbooks, where for real 
eigenstates one knows that 

i r+oo 

MP) = u n (x)e-^l h dx (24) 

j r+co 

= —j= I u n (x) [cos(px / h) — i sm(px / h)] dx (25) 
= A(p) - iB(p) (26) 

where A(p),B(p) are then purely real functions. This immediately implies that 

\M-P)\ 2 = \A(p)+iB(p)\ 2 = [A{p)f + [B{p)f = \A(p)-tB(p)\ 2 = \M+P)\ 2 - (27) 

The general expression for a wave packet solution constructed from such energy eigen- 
states is 

oo 

if>(x,t) =J2 a nMx)e~ iEnt/h (28) 

n=l 

where the expansion coefficients satisfy J2 n \a n \ 2 = 1. The time-dependent Wigner distribu- 
tion for this case is then 



1 P+OO 

P$\x,p;t) = —J r(x + y,t)ij(x-y,t)e 2ip y/ h dy 



OO OO r 

i(E m -E n )t/h ! 



n 6 



m=l n=l 

OO OO 



+oo 

irk 



u m {x + y) u n {x-y) e 2ipy/h dy 



= EEW^^^^'f^) (29) 

m=l n=l 

where, in general, we must calculate both diagonal (m = n) and off-diagonal (m ^ n) terms 
of the form 



P^' n \x,p) = / Um (x + y)u n (x-y)e 2i ^ h dy. (30) 



-oo 



While P$\x,p; t) is real, the off-diagonal Wigner terms are not necessarily so, but do satisfy 



P^ n \x,p) =P^ m >(x,p). (31) 



)(n,m) 
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III. INFINITE SQUARE WELL RESULTS IN POSITION AND MOMENTUM 
SPACE 



Because the Wigner function can be evaluated using either position- or momentum-space 
wavef unctions, we review the properties of, and interconnections between, these solutions 
for the infinite square well (ISW). The standard problem of a particle of mass m confined 
by the potential 

{0 for < x < L 
(32) 
oo otherwise 

has energy eigenvalues and position-space eigenf unctions (which are non-zero only within 
the well) which are given by 

h 2 ir 2 n 2 pi j2 /rmx\ 

n = IfrrfX = 2^ ( pn = n7[h ' L > Wlth Un ^ = V 1 sm VT~ ) ' ^ 

The position-space eigenf unctions have a generalized parity property about the midpoint of 
the well since 

u n (L -x) = ^1 sin ( n7r(L ~ X) ^j= -y|cos(n7r) sin (^) = (-1)" +1 u n {x) . (34) 

As with any such system, the eigenfunctions can be made orthonormal, which in this case 
can be readily checked explicitly by direct calculation of 



<«*K> = Il)J sm {—) sm {—) dx 

sin[(m — n)ir] sin[(m + n)7r] 



(m — n)n (m + n)n 

The momentum-space eigenfunctions are given by Eqn. (J2J) as 
1 



(35) 



»(p) = I" u n {x)e-^' h dx 

\J2ixh Jo 



e +in7r/2 sm[(pL/fi-mr)/2] _ ^_ in7v/2 sm[(pL/H + rar)/2] 



(pL/h — mr) (pL/h + nn) 



and the resulting probability density is given by 

/ T \ r„i„2l7„r /fc „_\ /ol „;„2 
|2 I L \ 



\Up)V = 



- 



sm 2 [(pL/h-nn)/2] + sin 2 [(pL/H + rm)/2) 



Kit ) \ {pL/h — rni) 2 (pL/h + nn) 2 

sm[(pL/h — mr)/2] sin[(pL/h + mr)/2] 



—2 cos(n7r) 



(pL/H — n7r)(pL/H + nir) 



(36) 



which will be useful for visualization purposes. 

While these forms demonstrate more explicitly the strong peaking of the probability 
amplitude near the expected values of p = ±p n = ±nnh/L, they do not make clear the finite 
extent (limited to the range [0, L]) of the position-space amplitude. In order to exemplify this 
dependence, we can evaluate the inverse Fourier transform, via Eqn. to explicitly obtain 
u n (x). Because such calculations will be useful in the evaluation of the Wigner distribution 
using momentum-space wavef unctions, we examine this seldom-discussed analysis |9] in some 
detail. We require 

i r+oo 

u n (x) = -== / Mp)e ipx/h dp 
\/2txTi J-oo 

L f c+ W2 f +0 ° c fefa-£/2Vft sm[(pL/h-n7r)/2] ^ ^ 



-i, 



2-K 2 h 2 { J _ oa (pL/h — nn) 



-oc 



_ e -inn/2 I e i P (x-L/2)/h sm[(pL/h + nir)/2} 
1 ' (pL/H + nir) 



oo 



" H)V^{ e+W2 ^- e - W2 ^}- 

Each integral, I a-, Ib, can be done in turn using a change of variables, giving, for example, 

M*-m/n sm[(pL/H-mr)/2) 



e wx-wn dp 
(pL/n — rm) 

— \ P + in7TX / L p-inn/2 [ + Sm(g) imq 



oo 



LJ J-oo q 

= g e +inirx/L e -in,/2 J °° ^) [ cog ( mg) + % sin ( mg )] dq 

where we have used 

1 fpL \ , (2x-L) 
q = — I — nn I and m = . (39) 

Z \ tl J Ij 

The piece of I a which includes the sin(mg) term vanishes for symmetry reasons (odd in- 
tegrand over a symmetric interval), while the component with cos(mg) can be done using 
standard handbook |45[ results, which we review in Appendix EI The result is 

ii for m 2 < 1 

" + °° sin(g) 



q 



cos(mg) dq 



tt/2 for m 2 = 1 (40) 
for m 2 > 1 



The restriction on m corresponds to 

m 2 = I =r_Jl ) < i _i < <+1 0<x<L (41) 



2x-L\ 2 2x-L 
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which explicitly demonstrates the finite extent of the position-space wavefunction which is 
non- vanishing only in the range [0, L], as expected. Performing the second integral (Ib) 
in the same way, and combining factors, we find that the non-vanishing position-space 
wavefunction is given by 



2n 2 h 2 V L 



^ ^g+in7r/2g+m7rx/L ^— inir/2 ^— inir/2^— inirx/L ^+inir/2 ^ 



\ sin (VH^.\ f or o < x < L (42) 



U n {X) 

L V L 

again, as expected. For future notational convenience, we can write this in the form 

u n {x) = ^1 sin K(x; 0, L) (43) 

where 



7Z(x; a, b) 



for x < a, x > b 

1/2 for x = a,x = b ■ (44) 

1 for a < x < b 

(We note that the results from Eqns. (}4T)j) . (}4*3*j) . and (pHj) can also be expressed in terms of 
of the Heaviside step-function, 0(£) ; as ; f° r example in Ref. 

Just as in position space, it is possible to explicitly demonstrate the ort honor mality of the 
momentum-space eigenfunctions, by calculating {4> m \4>n)- Since similar methods will also be 
useful in what follows, we illustrate one typical step in such an evaluation. One required 
integral, for example, is given by 

j = (Jl) e *(»-mWa [ +Q ° sin[(pL/ft - mn)/2] sm[(pL/h - nvr)/2] 
\7th J J_ 00 (pL/h — rrm){jpL/h — nir) 

The denominator can be written in a form which allows use of standard integrals, namely 



1 



1 



(pL/h — mn) (pL/h — nir) 



(46) 



(pL/h — mTc)(pL/h — rnr) (m — n)7t 

Appropriate changes of variables and integrals as in Eqn. (|40|) then give simple closed form 
expressions, and the complete result for (<f) m \<f) n ) is the indeed same as in Eqn. ([3*5]). 

Finally, for eventual comparison to the Wigner distributions for the ISW, we plot some 
standard representations of the position- and momentum-space probability densities for two 
low-lying states (n = 1, 10) in Fig. [TJ 
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IV. CLASSICAL PHASE-SPACE PICTURE OF THE INFINITE SQUARE WELL 



Phase-space plots of the motion of classical mechanical systems are increasingly stressed 
in undergraduate textbooks on the subject, perhaps because of their utility in the analysis 
or visualization of classically chaotic systems. The correlated time dependence of x(t) and 
p(t) in such systems is indeed easily and usefully visualized by parametric plots in the x — p 
plane, and the most familiar example is perhaps the case of the harmonic oscillator. For 
this case, the most general form of the solution can be written in the form 

x(t) = xa cos(ut + (p) and p{t) = —mux a sin(a;t + (p) = pa sin(a;t + (f>) (47) 

and the constant value of the total energy is 

isr +H^i 2 = B =^ = i- (48) 

This defines the elliptical phase-space path followed by the particle, and a classical descrip- 
tion of the joint probability densities requires that the particle be restricted to this path, 
namely, that 

Pcl{x,p) oc 6{E - p 2 /2m - mu V/2) . (49) 

For the classical infinite square well, the situation is at once conceptually simpler, but 
somewhat more subtle mathematically. The classical paths consist of simple 'back-and-forth' 



motion, at constant speed (p = ±p = ±\ / 2mE) with sudden (discontinuous) jumps at the 
walls (x = 0,L). This is illustrated in Fig. El with (solid) horizontal lines at ±po connected 
by (dashed) vertical 'jumps' at the walls. The classical joint probability distribution can 
then be described intuitively as 

P CL {x,p) oc 8{E -p 2 /2m) for < x < L (50) 

and the corresponding classical probability densities for position, Pcl{x), and momentum, 
Pcl{p), are obtained by integrating over one variable or the other, as in Eqns. (J7J) and 
(JHJ). For the p-space probability density, since there is no explicit x dependence inside the 
5-function, one finds that 

Pcl(p) oc / 5(E - p 2 /2m) dx oc 5{E - p 2 /2m) oc [6{p - p ) + 6{p + p )] . (51) 
J 2p 
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When properly normalized this reduces to 

Pcl{p) = g [Sip - Po) + <*(p + Po)] (52) 

or equal probabilities of finding the particle with p = ±po- The corresponding position-space 
probability 

Pcl(%) oc y 5{E — p 2 /2m)dp oc constant independent of x (53) 
for < x < L, or when properly normalized gives 

P CL (x) = - for < x < L . (54) 

Li 

These classical values are also shown in Fig. El (dashed horizontal line for Pcl{x) and two 
vertical dashed 'spikes' for Pcl(p)) to be compared to a large quantum number (n = 10) 
solution of the quantum case for comparison. 

V. WIGNER DISTRIBUTION FOR THE INFINITE WELL: EIGENSTATES 

For the calculation of the Wigner distribution for energy eigenstates in the ISW, we will 
first work in position space and use 

p w(x>P) = ^ / u n{x + y) u n {x - y) e 2 ^ h dy (55) 

since the position-space eigenstates, u n (x), can be made purely real. The limits of integration 
are determined by the restriction that the u n (x ± y) are non- vanishing only in the range 
[0, L], so that we must simultaneously satisfy the requirements 

0<x + y<L and < x - y < L . (56) 

This leads to upper and lower bounds for the integral over y in Eqn. (|55|) which depend on 
x via 

- x < y < +x for < x < L/2 (57) 

-(L-x) <y <+{L-x) for L/2 < x < L. (58) 
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Thus, over the left-half of the allowed x interval, [0, L/2], we have 



w 1 



x,p) 



1 

irh j _ ,. 
2 



sm 



nn(x + y) 
L 



sm 



n7r(x — y) 
L 



2ipy/h 



dy 



sm[2(p/h — nir/L)x] sin[2(p/h + nir/L)x] 



4(p / h - n-K j L) 4(p / h + nit / L) 

( 2nirx\ sm(2px/H) 1 
^"^J (2p//i) J 



(59) 



— cos 



while over the right-half of the interval, [L/2, L], one makes the replacement x — > L — x. This 
form has been derived before [^| , Q| , [2^ although in at least one reference it is written in 
terms of Bessel functions {jo(z)) which somewhat obscures its simple derivation. One can 
then demonstrate explicitly that 

"L r+00 

P^\x,p)dp=\u n {x)\ 2 (60) 



P$\x,p) dx = |0 n (p)P 



and 



using standard integrals, or ones described in Appendix |XJ 

The evaluation of the Wigner function using momentum-space wavefunctions, using 
Eqn. (J5J) integrated over all p-values, is also instructive, as it naturally leads to the same 
restrictions on x as seen in Eqn. (J41|) . as well as the x — > L — x identification over the two 
half-intervals, as in Eqn. (|58jl . We exhibit some of the relevant parts of this calculation in 
Appendix El 

With these results in hand, we exhibit examples of the Wigner function for the ISW, for 
the n — 1 ground state (Figs. ElandEJ) and for an excited (n = 10) state (Figs. El and EJ), 
using h = L = 1 for simplicity. For ease of visualization, for both cases we show separate 
plots for Pyy(x : p) > and P^\x,p) < values. For the ground state, the projection of 
Pw\ x ->P) on ^° ^ ne x ~ an d p-axes to obtain |w n (a;)| 2 and \<p n {p)\ 2 (as in Fig. EJ via Eqn. (f6T)j) 
is straightforward to visualize, and we also note that this relatively smooth ground state 
wavefunction gives a Wigner distribution which is almost everywhere positive (Fig. OJ), but 
with a small negative contribution as seen most clearly in Fig. |3J 

For the n — 10 (or any n » 1) case, the structures are more surprising, but still yield 
the appropriate probability densities upon projection onto the x- and p-axes. The obvious 
'triangular' form of both the 'fin'-shaped features along the p = ±p n axes and the central 
'spines' along the p = axes can be easily seen to arise from the form in Eqn. (|59jl in those 
limits. For those cases, the sin [2 Fx/ ft] /F form (where F = (p — p n ), (p + p n ), or p) of each 
term (for < x < L/2) gives a linear dependence on x when F — > 0, while for L/2 < x < L 
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this is replaced by the L — x factor, giving the triangular dependence obvious from Fig. |SJ 
The smooth features along the p = ±p n axes (seen only for positive values of P^\x,p)) 
are obviously reminiscent of the classical phase-space picture in Fig. |2l while the highly 
oscillatory structures 'inside' the classical rectangular boundary (obvious in both Figs. 
andEJ) locally average to the small, but non- vanishing values of the cross-term in Eqn. (J36|) 
and are just as clearly purely quantum mechanical in origin. 



VI. WIGNER DISTRIBUTION FOR THE INFINITE WELL: WAVE PACKET 
SOLUTIONS 



For a general wave packet solution for the infinite square well, we require the on- and 
off-diagonal terms in Eqn. Using the position-space eigenstates in Eqn. (JBH1), over the 

interval [0, L/2] we find that 



P 



(m,n) 



W 



x,p) = — 
irn 

1 



rrnr(x + y) 
L 



mr(x 



y) 



L 



-,2ipy/h 



dy 



3 +i(m— n)irx/ 



L sm[(2p/h + (m + n)ir/L)x] 



+e 



(61) 



2pL/h + (m + n)ir 
-i{m-ny x /L sm[(2p/n - {m + n)n / L)x] 
2pL/h — (m + n)7r 
_ e+i{m+ n)nx/L sin[(2p/h+ (m-n)n/L)x} 

2pL/h + (m — n)n 
_ e - i{m+ n)nx/L sm[(2p/h-(m-ri)7r/L)x\ 
2pL/h — (m — n)ir 

and it is easy to check that this result reduces to the expression in Eqn. (|59|) when m = n. 
In order to extend this to the interval [L/2, L], it is important to note that the substitution 
x — > L — x should be made only in those terms arising from the integration over dy, namely, 
the sm[(2p/h± (m ± n)/L)x] terms. 

These results can then be used to evaluate the time-dependent Wigner function for any 
initial wave packet in the infinite square well, with the simplest non-trivial example being 
that for a two-state system, ip(x, 0) = 7«i(x) + 5v,2(x), which can then be compared to 

|. Such states, 



more standard representations in position space or momentum space 
while exhibiting non-trivial dynamical behavior (sometimes said to mimic time- dependent 
radiating systems 471]). are still highly quantum mechanical and to approach the semi- 
classical limit, we require more localized states, constructed from higher quantum number 
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eigenfunctions. 

We can examine the time dependence of such an initially localized state by choosing a 
standard Gaussian of the form 

ljj {G) ( X , 0) = —L= e -(x-*o) 2 /2bl e ip (x-x )/h ^ 



where we will always assume that xq is chosen such that ip(G)(x, 0) is sufficiently contained 
within the well so that we make an exponentially small error by neglecting any overlap with 
the region outside the walls, and may thus also ignore any problems associated with possible 
discontinuities at the wall. In practice, this condition only requires the wave packet to be a 
few times Axo = away from an infinite wall boundary. With these assumptions, we 

can then extend the integration region from the finite [0, L] interval to the entire 1-D space, 
giving the (exponentially good) approximation for the expansion coefficients 4sj], j3] 

L 

[u n (x)} [iI>g(x, 0)] dx 

+ 0O 

[u n {x)\ [iI>g(x, 0)] dx (63) 



— | / j^H I innxp/L -b 2 (p +mrh/L) 2 /2h 2 _ -innx /L -b 2 (p -mrh/ L) 2 /2h 2 l 

The position-space and momentum-space wavefunctions are then given by 

oo oo 

i/>(x, t) = Y, a nu n {x) e-^l h and 0(p, t) = ^ a n <p n (p) e^"^ (64) 

n=l n=l 

while the time-dependent Wigner distribution is given by Eqn. (|29J). 

The time dependence of any general quantum state (not necessarily the ISW) is deter- 
mined by the exp(—iE n t/h) factors, and for highly localized, semi-classical wave packets, 
where one typically can expand E(n) about a (large) central value of the quantum number, 
n » 1, one can write this dependence in the form 4l|, 42 1 

e -iE n t/n = exp ^_ z /ft E(n )t+ (n-no)E'(n )t+ ^(n - n ) 2 E" (n )t 

+ h n -n fE'"(n )t + -- 
6 

= exp (— iu> t — 2ni{n — n )t/T c [ — 2-aiin — n ) 2 t/T rev (65) 
-2m(n - n ) 3 t/T super H ) 
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where each term in the expansion (after the first) defines an important characteristic time 
scale, via 

2irh 2ixh 2ixh 

d ~ \E'{n )y rev ~ \E"{n )\/2' ^'W^U^ 1 

namely, the classical period (T c i), the revival time (T rev ), and the super-revival time (T super ). 
The classical periodicity for the ISW in this formalism is given by 

2irh 2mL 2 2L 2L , p no , . 

Td = 1^77 r, = t = 777 /rw 1 = — where v no = — 67 

\h'(no)\ fmriQ [(mrno/ L)/m\ v no m 

and v no is the analog of the classical speed, giving a result in agreement with classical 
expectations. Depending on the relative values of T c i and the spreading time, to, there can 
be many obvious similarities to the classical motion, easily visualized through an expectation 
value analysis j^]. 

As an example of quasi-classical time evolution of a quantum state, we consider a Gaussian 
wave packet characterized by the parameters 

2m = L = h = 1 A.x = 4= = — « 1 = I* (68) 

y/2 20 1 ; 

= = 40vr » io = A Po = -J- x = L/2 = 0.5 . (69) 

L 2Axq 

We use a central values of Uq = 40 so that the ratio of classical periodicity to spreading time 
is T c i/to — 10/-7T and significant spreading can be seen even over half a classical period. 

We first show in Fig. [7| the position- and momentum-space probability densities for the 
initial wave packet (t = 0, dashed) and half a classical period later (t = T c ;/2, solid). 
The spreading in position space is clearly visible, while the switch from momentum values 
centered around +po to — po after the first 'collision' with the infinite wall is also apparent. 
For comparison, we show Pyy(x,p;t) for the same two times in Fig. |H] and note the two 
highly localized 'lumps' of probability, centered at the appropriate locations in phase space, 
but with the obvious spreading in position space, consistent with the form in Eqn. (|17jl. 
and with positive-definite values in each case. This behavior is similar to the classical phase 
space picture in Fig. El for the two states denoted by squares, separated by T c i/2. 

For comparison to this quasi-classical behavior, we also show results for t = T c i/4, corre- 
sponding to the time when the classical particle would 'hit' the wall, both for the position- 
and momentum-space probability densities (in Fig. ^ as well as for the Wigner function 
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(in Fig. ITUll . In this case, the description of the time-dependent solution in terms of sums 

n n n 

of image states |48j . j51lj . |52j is appropriate, with the large, interference term at p ~ in 
the Wigner distribution arising from effects such as in Eqn. (|19|) . As discussed above, the 
Wigner function is positive-definite only at times when it can be well approximated as a 
single isolated Gaussian, such as at t = and t = T c i/2 in Fig. |Sl or, as we will see, at 
the revival time T rev . As seen here at T c //4, at many other times the Wigner function can 
become negative and we will henceforth only show the Pw(x,p;t) > values for ease of 
visualization. 

For longer time scales, we require the revival time, T rev , which is given by 

T rev = , 2 ** u = ™ = ^ = (2n )T d (70) 
| E" (n ) |/2 E hi: V 0J y 1 

and for the infinite square well no longer time scales are present due to the purely quadratic 

dependence of the energy eigenvalues. The revival time scale can clearly be much larger 

than the classical period for wave packets characterized by no >> 1 (giving T rev /T c i = 80 

for the examples presented here.) 

The quantum revivals for this system are exact since we have 

+oo +00 

^(s, t + T rev ) = a n u n (x)e- iE ^ t+T ^ h = £ a n u n (x)e-^/V i2 ™ 2 = ^(z, t) (71) 

n=l n=l 

and any wave packet returns to its initial state after a time T rev . At half this time, t = T rev /2, 
the wave packet also reforms 53^ since 

oo 

ip(L -x,t + T rev /2) = a nUn{L - x ) e - iE " t / h e - iE " T -/ 2h 

n=l 

oo 

= ^a n u n {x)e- lE - t/h \{-l) n+l e- in2 A (72) 

n=l 

where we have used the symmetry properties of the u n (x) from Eqn. (|34|1. This implies that 

\iP(x,t + T rev /2)\ 2 = \i;(L-x,t)\ 2 (73) 

so that at half a revival time later, any initial wave packet will reform itself (same shape, 
width, etc.), but at a location mirrored about the center of the well. Using Eqn. (|72|). one 
can also show that 

|0(p,t + T_/2)| 2 = |</)(-p,t)| 2 (74) 
18 



so that the packet is moving with 'mirror' (opposite) momentum values as well, and therefore 
in the opposite corner of phase space. 

More interestingly, at various fractional multiples of the revival time, pT rev /q, the wave 
packet can also reform as several small copies (sometimes called 'mini-packets' or 'clones') of 
the original wave packet, with well-defined phase relationships. (The original mathematical 
arguments showing how this behavior arises in general wave packet solutions were made by 



Averbukh and Perelman in Ref. 



). For example, near t ~ T rev /4 the wave packet can be 



written linear combination of the form 

^(x, t » T rev /A) = ±= [e~™ /4 ^ cl (x, t) + e +M/4 ^ cl (x, t + T d /2)] 

where tp c i(x,t) describes the short-term time development of the initial wave packet ob- 
tained by keeping only the terms linear in n in the expansion in Eqn. (|65|). The position- 
and momentum-space probability densities at T rev /4 are shown in Fig. ^2 illustrating this 
phenomena, but this behavior is much more interestingly visualized through the Wigner 
function description in Fig. El The Wigner function in this case is well-represented by 
the terms in Eqn. (|19|h where the cross-term is only oscillatory in the p variable, since the 
corresponding values of xa = %b = L/2 = xo are the same in this case. 

A similar, but even richer, situation can be seen at T rev /3, where there are three dis- 
tinct 'mini-packets' in position space, and corresponding interesting interference structure 



in momentum space, as shown in 
fractional revival time is given by 



ig- El The explicit form of the wavefunction near this 



il>(x, t » T rev /3) = ty d ( Xj t ) + e 2m/s {i() d {x, t + T c/ /3) + ij> d (x, t + 2T d /3)}] (75) 

and exactly this type of highly correlated state is obvious in the Wigner function visualization 
in Fig. El with three smooth 'lumps' corresponding to the diagonal terms in Eqn. (|19p. and 
three oscillatory 'cross-terms' at the average value locations in phase space predicted by 
Eqn. (J2Uj) . with the appropriate 'wiggliness'. 

Finally, at other 'random' times during the time evolution of such states, not near any 
resolvable fractional revivals, the wavefunction collapses to something like an incoherent sum 
of eigenstates, with little or no obvious correlations. A view of this behavior at such a time 
(represented here by t = T* = 16T rev /37) using \ip(x,t)\ 2 and \4>(p, t)\ 2 is shown in Fig. ITol 
as well as using the Wigner function in Fig. El 
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VII. CONCLUSIONS AND DISCUSSIONS 



While the infinite well is one of the most standard model problems in all of introduc- 
tory quantum mechanics, it continues to be used as a benchmark system for the analy- 
sis and visualization of new quantum phenomena, including wave packet revivals. In a 
similar vein, the Wigner function, which was invented more than 70 years ago, continues 
to play an important role in the analysis, understanding, and visualization of quantum 
systems, and in related fields such as signal analysis and quantum optics. We have pro- 
vided a thorough review of both the analytic structure of the Wigner function for energy 
eigenstates as well as for time-dependent Gaussian wave packet states for this important 
exemplary quantum system, focusing on the visualization of wave packet states at frac- 
tional revival times. (Additional images not included in this paper can be accessed at 
http : //webphysics . davidson . edu/mjb/wigner/ ) 



Given the relative straightforwardness of the calculations involved here for the ISW, it is 
easy to imagine extending the results presented here to simple extensions, such as a finite 
square well, or an asymmetric infinite well of the form 



V{x) 



+00 for x < — b 

for -b < x < 

(76) 

+V for < x < +a 
+00 for +a < x 



This last form is interesting as the eigenstates exhibit less trivial correlations between the 
magnitude/wiggliness of the position-space wavefunction (due to the varying classical speeds 
in the two sides of the well) which can then be connected to the behavior in momentum space 
more directly using the Wigner function. This example is also instructive since, because of 
the unphysical discontinuity of the potential, there are typically surprising results [5^ when 
one compares quantum results to classical expectations for position- and momentum-space 
probability densities. This type of asymmetric well is also of the form for which coherent 
charge oscillations have been observed experimentally consisting of a two state system, 

V(x, t) = 7 [u n (x)e-^ h ] + 6 [u b (x)e- lE ^ h ] (77) 
with E a < V and E b > V . 
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APPENDIX A: TRIGONOMETRIC INTEGRALS 



For many of the calculations included here, we require versions of the integral 



sm(z) cos mz 



dz = < 



for m < — 1 and m > +1 
7r/2 for m = ±1 
7i for m 2 < 1 



(Al) 



which is a handbook |45( result. We can briefly justify these results, starting with the single 
integral 



+oo 



sin z 



dz = 71 



(A2) 



which itself can be derived using contour integration. This can be generalized to 



sin (mz\ 



dz = < 



+7T for m > 
for m = 
—n for m < 



(A3) 



by a change of variables, and considering the special m = case separately. The general 
integral in Eqn. (jAl|) can then be obtained by writing 

1 



sin(z) cos(mz) = - {sin[(l + m)z] + sin[(l — m)z]} 

and using Eqn. (jA3|) twice. The special case of m = 1 is done by noting that 

^ +00 sin(2z) . 1 r+co 



(A4) 



sin m cosu 



dz 



dx — — 
2z 2 



svayw) , 7i 

dw = — 

w 2 



(A5) 



where w = 2z and a simple trig identity is used. Another integral that is useful for normal- 
ization calculations is 



+°° sm 2 (z) 



dz = 71 



(A6) 



which is another handbook result derivable using contour integration techniques. 
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APPENDIX B: WIGNER DISTRIBUTION FROM MOMENTUM-SPACE 
WAVEFUNCTIONS 



The evaluation of the Wigner function, using momentum-space wavefunctions as in 
Eqn. ©, requires the integral 



+00 



dq e 



-2iqx/h +i(p+q)L/2h -i{p-q)L/2h 



(Bl) 



x < e 



x < e 



-»W2 sin [((P + q)L/h-nir)/2] _ c+in7r/2 sin[((p + q)L/h + nn)/2] \ 

[(p + q)L/h - nn] [(p + q)L/h + nn] J 

-inn/2 sm[(ip-q)L/h-rnr)/ 2] _ c _ iW7r/2 sin[((p - q)L/h + nn)/2] 1 

\{p — q)L/h — nn] [(p — q)L/H + nn] J 

and we briefly sketch out some of the necessary steps in the evaluation of P$ p) i n this 

approach. 

Combining the various complex exponentials and using the fact that the Wigner function 
must be real, we are left with integrals such as 



cos 



qL 

h 



2x 



sin[((p + q)L/h — mr)/2] sin[((p — q)L/H — nn)/2] 



[(p + q)L/H— n-rc] 



[ip — q)L/h — nn] 



In this case, the appropriate partial fraction identity to rewrite the denominators is 



1 



((P + q)L/h — nn)((p — q)L/h — nn) 2(pL/h — nn) { ((p + q)L/H — nn 



1 



dq . 
(B2) 

(B3) 



((p — q)L/K — nn) 

and the remaining integrals can be done using variations on Eqn. (|40p. Upon combining 
various factors, one obtains integrals of that form not only with m = (L — 2x)/L, as in 
Eqn. (JSHJ), but similar ones with m = (L — 4x)/L and m = (3L — Ax)/L. These terms give 
rise to limits on the x dependence of the form JZ(x; 0, L/2) and 7Z(x; L/2, L) where TZ(x; a, b) 
is defined in Eqn. (144)1 . For example, one intermediate result can be written in the form 



T 



cos 



— sin 



pL 



nn 



pL 
T 



nn 



L-2x 
L 

L - 2x 







sin 





cos 



pL 



nn 



n(x;0,L) (B4) 
{TZ(x; 0, L/2) - TZ(x; L/2, L)} 



which, in turn, gives 



r 







for x < 0, x > L 
sin[(2(p/h — nn/L)x] for < x < L 

sm[{2(p/h - nn/L){L - x)\ for L/2 < x < L 



(B5) 
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so that the 'split' definition of (x,p) in the two half intervals arises very naturally and 
the complete result in Eqn. is reproduced using momentum-space methods, including 
the non-trivial x dependence. 



APPENDIX C: PROBLEMS 



PI: Complete the proof that the <f) n (p) for the ISW are orthonormal by explicit calculation 
of {(f) m \4>n), completing the steps in Sec. and making use of identities such as Eqn. 

P2: Using results from any quantum mechanics textbook, evaluate the Wigner distribution 
for the ground state and first excited state of the simple harmonic oscillator. Answer: The 
harmonic oscillator eigenstates can be written in the form 

u n (z) = A n H n (z)e~ z2/2 where A n = 1 and z = ^ = 7= (CI) 

and the H n (z) are the Hermite polynomials of order n. The Wigner functions for n = 0, 1 
are then given by 

1 



P$\x,p) = —e-o 2 (C2) 
7rn 

P&\x,p) = \^ 2 (2p 2 -l) (C3) 



where 



This is one of the simpler results which shows explicitly that the Wigner function need not 
be positive-definite. Note that the first of these results is consistent with the expression for 
the free-particle Gaussian in Eqn. (JT7j) for vanishing xq, po, and t. The Wigner function for 
arbitrary n has been evaluated (see, e.g., Ref. jsj) with the result that 

P { w ] = ^V 2 L n (2p 2 ) (C5) 
where L n (z) are the Laguerre polynomials. 

P3: Show how the phase-space plot for the 1-D harmonic oscillator can be used to generate 
the classical probability densities for position- and momentum-variables, namely, show how 
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one obtains 

P Ci (x) = — ^= and P CL (p) = —7^^ ■ (C6) 

Partial answer: The classical probability distribution from Eqn. (J49|) can be written in 
the form 

P CL (x,p) oc 6(p(x) 2 -p 2 ) = ^-^[6(p(x)-p) + 6(p(x)+p)} (C7) 

where 



p(x) = \JlmE — m 2 uj 2 x 2 = muy x\ — x 2 . (C8) 
Integration of Eqn. (jC7|) over dp then gives 

PcL<x) K W) K (C9) 

which when properly normalized (integrated over the classically allowed interval (—xa, +%a)) 
gives the first term in Eqn. (|C6I) . 



P4: Using either the momentum- or position-space wavefunctions in Eqns. f!13|) or (fT4*|) for 
the free-particle Gaussian wave packet, show that the Wigner function is of the form in 
Eqn. (JT71). 

P5: Complete the proof that the result for the Wigner function for eigenstates of the ISW 
in Eqn. can be obtained using momentum-space wavefunctions, using the techniques in 
Appendix |Bj 

P6: Use the results in Eqn. f!29|) and (|62|) to write down the time-dependent Wigner function 
for the simple two-state system in the infinite well 

il>(x,0) = — [ Ul {x)+u 2 {x)) (CIO) 
v2 

and generate plots of (x, p; t) for various times. Compare the results to standard images 
of the position-space and momentum-space probability densities for this problem 
What is the only time scale associated with this system? and does it have anything to do 
with a classical periodicity? 
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Iu n (x)l 2 vs. x I4> n (p)l 2 vs. p 




-Pio + Pio 
FIG. 1: Plots of the position-space probability density, |-u n (x)| 2 versus x, (left) and momentum- 
space probability density, \4> n {p)\ 2 versus p, (right) for energy eigenstates in the infinite square well 
for n = 1 (top) and n = 10 (bottom) cases. The horizontal dashed lines on the left correspond to 
the (classical) flat probability density given by Pcl(x) = 1/L from Eqn. (|54|). while the vertical 
dotted lines on the right correspond to the 5- function classical distribution in Eqn. (|52|l. 
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x-p plane 




FIG. 2: Classical phase-space picture of solutions of the infinite square well. The classical phase- 
space trajectory has the particle moving with constant speed (momenta given by ±po) between 
the walls at x = 0, L (solid horizontal lines) with discontinuous changes in velocity (momentum) 
due to the collisions with the walls (dashed vertical lines). Projections onto the x- and p-axes give 
the classical probability densities in Eqns. and l[54|). These are compared with the quantum 
counterparts, Pqm{x) = \u n (x)\ 2 versus x, and Pqm{p) = |0n(p)| 2 versus p, for the case n = 10. 
The pairs of points in phase space indicated by squares (stars, diamonds) are separated in time by 
half a classical period. 
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FIG. 3: Plot of the Wigner distribution from Eqn. H59j) for the n = 1 energy eigenstate in the 
infinite square well. Only the positive (pff(x,p) > 0) parts are shown. 



32 




I 



FIG. 4: Same as Fig- El but with only the negative (— Pyy (%,p) > 0) parts shown. This shows 
that the Wigner function for the ground state of the infinite well is almost, but not quite, positive 
definite. 
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I\|/(x t)l 2 VS. X |(|)(p t)l 2 vs. p 




L -P +P 

FIG. 7: Plots of the position-space probability density, |^(x,£)| 2 versus x, and momentum-space 
probability density, \(f>(p, t)\ 2 versus p, for a Gaussian wave packet solution in the ISW. Times 
corresponding to t = (dashed) and T c i/2 (solid) are shown. The parameters of Eqn. (|69|) and 
no = 40 are used. 
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FIG. 8: Plot of the Wigner function, F\y(x,p;t) versus (x,p) as a function of time for t = and 
t = T c i/2, to be compared to Fig. [7| 
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I\|/(x t)l 2 VS. X 



l(|)(p t)l 2 vs. p 
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0.5L L -P +P 

FIG. 9: Same as Fig. [7J except for t = (dashed) and t = T c i/4 where the classical particle 
would be hitting the wall. (Recall that the two momentum peaks for the 'collision' time are not 
symmetrically placed at ±po |53] since the high-momentum components of the initial wave packet 
arrive at the wall, and hence are also reflected, first.) Note that the \ip(x,t)\ 2 is plotted over the 
interval [L/2, L] to show the collision with the wall in more detail. 
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0,5 

FIG. 10: Same as Fig. |H1 but for t = T c i/A where the classical particle would be hitting the wall. 
Only positive values (P\y(x,p;t) > 0) are shown and the x interval [L/2,L] is used, as in Fig. El 
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I\|/(x t)l 2 VS. X 
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FIG. 11: Same as Fig. [7[ except for t 
(dashed) 



L -P +P 

(dashed) and for a fractional revival at t = T rt 
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FIG. 12: Same as Fig. El but for the fractional revival at t = T rev /4, to be compared to Fig. 1111 
Only positive values (Pw(x,p;t) > 0) are shown. 
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FIG. 13: Same as Fig. [7J except for t = (dashed) and for a fractional revival at t = T rev /3 (solid). 
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FIG. 14: Same as Fig. |H1 but for the fractional revival at t = T rev /3, to be compared to Fig. EH 
Only positive values (Pw(x,p;t) > 0) are shown. 
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FIG. 15: Same as Fig. [7J except for t = (dashed) and for a general time, T* = 16T rev /37 (solid), 
during the collapsed phase, not near any resolvable fractional revival. 
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FIG. 16: Same as Fig. IH1 but for a general time (T* = 16T rev /37) during the collapsed phase, not 
near any resolvable fractional revival, to be compared to Fig. El Only positive values (Pw(x, p; t) > 
0) are shown. 
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